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Abstract 

We derive a nonperturbative result for the Higgs propagator by calculating 
the NLO correction in the 1/N expansion of 0{N) sigma models. For dealing 
with the expansion beyond leading order, we develop techniques to cal- 
culate certain infinite sets of multiloop Feynman diagrams. We then calculate 
the Higgs lineshape in fermion scattering, as seen for instance at a muon col- 
lider. We compare with the existing two-loop perturbative results. There is a 
considerable discrepancy between perturbation theory and the 1/N expansion 
in LO. By including the NLO correction, this discrepancy is reduced dramati- 
cally. The results are in very good agreement with the two-loop result even for 
masses as high as 850 GeV. A maximum of the mass is reached at 930 GeV. 
For phenomenological purposes, we give a simple approximate relation between 
the Higgs mass and width. 
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Abstract 

We derive a nonperturbative result for the Higgs propagator by calculating 
the NLO correction in the expansion of 0{N) sigma models. For dealing 
with the 1/N expansion beyond leading order, we develop techniques to cal- 
culate certain infinite sets of multiloop Feynman diagrams. We then calculate 
the Higgs lineshape in fermion scattering, as seen for instance at a muon col- 
lider. We compare with the existing two-loop perturbative results. There is a 
considerable discrepancy between perturbation theory and the 1/N expansion 
in LO. By including the NLO correction, this discrepancy is reduced dramati- 
cally. The results are in very good agreement with the two-loop result even for 
masses as high as 850 GeV. A maximum of the mass is reached at 930 GeV. 
For phenomenological purposes, we give a simple approximate relation between 
the Higgs mass and width. 

1 Introduction 

With the discovery of the vector bosons and of the top quark behind us, the next major 
question in the phenomenology of the weak interactions is the search for the Higgs 
boson. Within the standard model the properties of the Higgs boson are determined 
when its mass is fixed, as the tree level selfcouplings are proportional to rrijj. As 
long as the Higgs particle is light, there are no fundamental problems in determining 
its properties from perturbation theory. However, when the Higgs particle becomes 
heavy, 0{TeV), the selfcouphng becomes large, and perturbation theory becomes 
unreliable. 

Therefore one would like to find an approximation to Higgs physics beyond pertur- 
bation theory which is applicable for large couplings as well. When one realizes that 
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the Higgs sector of the standard model is nothing but an 0(4) hnear sigma model, 
the expansion in of the 0(A^)-symmetric sigma model suggests itself. 



The expansion has a long history — [11], and has also been applied to the 
Higgs boson of the standard model [jl], 3- The results in leading order in 1/iV are 
interesting. They indicate a saturation of the Higgs boson mass when the coupling 
increases. Also the relation between the Higgs width and mass is very different at 
larger couplings. However, the situation at leading order is not satisfactory because 
the lowest order is numerically very far from perturbation theory for small Higgs mass. 
Perturbation theory has recently even been extended to the two-loop level jT^ , p!3|, p^, 
and should be well under control for small Higgs mass. This shed for a long time 
doubts on the relevance of the 1/A^ approach. In order to clarify the situation it is 
therefore clearly desirable to know the next-to-leading order contribution in 1/A^. 

What makes the 1/A^ expansion particularly interesting is that it has certain fea- 
tures which are absent in perturbation theory, and which are to be expected from an 
exact solution. One feature is the absence in physical results of any renormalization 
scheme dependence. As we shall see, the next-to-leading order solution which we 
will derive has another such feature, namely the finiteness of wave function renormal- 
ization constants. Interestingly enough, this occurs after the summation of Feynman 
diagrams which are all ultraviolet divergent. 

The 0{N) sigma model at leading order in the 1/A^ expansion was discussed in 
detail by a number of authors. These works study the effective potential and the 
question of the occurrence of spontaneous symmetry breaking in 0{N) models, as 
well as two- and four-point Green functions. These aspects were studied from two 
points of view. One attempts to treat the 0{N) sigma model as a fundamental, 
renormalized theory; and the other as an effective theory, with a built-in cutoff. 
Although a large amount of work was done, most of it is limited to the leading order. 
Root PI wrote down expressions for the effective potential in higher order, but did not 
evaluate them explicitly. The reason for this is twofold. First of all, the expressions 
are quite complicated. More fundamental is the occurrence of a tachyon in the theory. 

A striking feature of the leading order solution is the presence in the propagators 
of a tachyonic pole. As long as one is concerned only with the leading order, and 
the Higgs coupling is not very strong, the causality violating effects induced by the 
tachyon are at least numerically negligible. They are suppressed essentially by the 
Landau scale. However, if one wants to calculate higher order corrections, these effects 
are unavoidable because the tachyon appears then in loops and leads to pathological 
solutions. For instance, it was shown by Root p| that the effective potential becomes 
complex in next-to-leading order. 

The tachyon problem is similar to the Landau ghost in perturbation theory. How- 
ever, the occurrence and position of the Landau ghost relies on using a perturbative 
expression precisely where one knows this expression is not a good approximation 
anymore. For the case of the tachyon problem in the 1/A^ expansion this is no more 
the case. The propagators in the 1/A^ expansion are summed to all orders in the 
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coupling constant. Thus their vahdity does not depend on the strength of the cou- 
phng. A possible interpretation for the presence of the tachyon is that it indicates 
the triviality of the theory. 

One simple and consistent way to deal with the tachyon is to regard the 0{N) 
model as an effective theory and to introduce explicitly a cutoff under the tachyon 
scale 0. However, we prefer to use in this paper a different approach. We introduce 
a scheme to subtract minimally the tachyon from the Green functions order by order 
in the 1/N expansion. This procedure is well-adapted to the calculation of higher 
order effects in 

The outline of this paper is the following. In section 2 we introduce the model 
and the leading order for fixing the notations. In the following section we discuss 
and motivate in detail the prescription for the treatment of the tachyon problem. 
In section 4 we describe the next-to-leading order calculation and the methods for 
evaluating the multiloop diagrams involved. In section 5 we treat the Higgs lineshape 
nonperturbatively and discuss the results. Section 6 summarizes the conclusions of 
the paper. 



2 Renormalization and leading order solution 

The 0{N) sigma model at leading order in the expansion was already discussed 
extensively in the literature. For this reason, we will not repeat here this discussion. 
In the following, for fixing the notations, we will only derive the two-point functions, 
assuming that the theory has a ground state with spontaneously broken symmetry. 
The starting point of the calculation is the following Lagrangian: 

c = ^a.$oa^$o - f - (1) 

Here $ is a scalar field with components i = 1, N. In principle one 
could use this Lagrangian in order to calculate contributions in different orders of 
1/N. However, the combinatorics becomes very complicated and it is advantageous 
to introduce an auxiliary field X) and add a nondynamical term to the Lagrangian 



according to a combinatorial trick proposed in ref. [11 1: 



This form of the Lagrangian has the same physical content as before. However, 
the Feynman rules are changed, thereby facilitating the counting of powers of 1/A^. 
For doing higher order calculations, the introduction of the extra field x is practically 
the only way not to get lost in the combinatorics - see also ref. . 
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The fields and coupling constants in eqns. 1, 2 are bare quantities. In order to 
perform renormalization, we make the following substitutions: 





(3) 



= -y(l + AO 



(4) 



X° = xZx + X + ^t^ 



l,...,iV-l 



(5) 
(6) 
(7) 



Here \/Nv and x the expectation values of the fields $ and^ in the broken 



symmetry ground state. For the case of the standard model Higgs sector, = 4 and 
V = 123 GeV. At^ and At^ are the tadpole counterterms corresponding to the fields 
a and 

The tadpole counterterms Ata- and At^ can be determined in the following way. 
In the original Lagrangian of eq. 1 only one tadpole counterterm is needed. Upon 
introduction of the auxiliary field x two such counterterms are present, but they are 
related to each other. Some contributions can be shifted from the definition of one 
counterterm to the other. It is then convenient for performing a next-to-leading 
order calculation to define the tadpole counterterms so that the vacuum expectation 
values of the fields a and x do not receive corrections in higher orders. By solving the 
gap equation in leading order, it is well-known that one finds in the broken symmetry 
ground state x = 0. By requesting that this simple relation be preserved in higher 
orders, we fix uniquely the relation between the two tadpole counterterms At^ and 
At^. Therefore, we define At^r and At-^ by the condition that the one-point functions 
for the fields a and x vanish. It is easy to convince oneself that this condition ensures 
automatically that the Goldstone theorem is valid, and the vr fields remain massless. 

Within perturbation theory, the counterterms defined above are power series of 
the renormalized coupling constant A. In the expansion they are power series in 
As it turns out, the contributions to most of these counterterms vanish in the 
leading order of the 1/N expansion: 
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i{N — l)a^^\s) = ( \ + counterterm 



Figure 1: The leading order bubble diagram. 

At, = ^*. + o(^ 

= ^«. + o(-^) (8) 

The only counterterm which is needed at leading order in is actually the ^A'^^-' 
counterterm of the quartic interaction. 

In order to calculate the propagators, one has to evaluate the era, xo", XX 
TTTT self-energy diagrams, and then to invert the matrix. As it turns out, the only 
nontrivial contribution to the proper self-energy functions of the theory in leading 
order of the expansion is due to the diagram shown in fig. 1. The ultraviolet 
divergence of the bubble integral is absorbed in the coupling constant renormalization, 
5\^^\ One finds the following well-known leading order propagators: 



D. 



s — m?(s) 

1 isni^{s) 
Nv"^ s — m'^{s) 
1 im?is) 



Nv s - m?'{s) 



1 

D^.^j{s) = 6ij- (9) 
s 

Here, the quantity m^(s) plays the role of an energy dependent mass, and is given 
by the bubble diagram of fig. 1 : 

where /i is the subtraction scale of the bubble diagram. 

We would like to reemphasize that, although a renormalization scale appears in 
eqns. 9, the physical predictions of this solution are free of any residual renormaliza- 
tion scheme dependence. This is because perturbation theory was summed up in all 
orders. Whenever one uses this solution of the 0{N) sigma model to relate physical 
observable quantities, the renormalization scale /i can be eliminated from the result. 
As an example, one can convince oneself that the relation between the Higgs mass 
and width is independent of the intermediary renormalization scheme used M. 
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3 Tachyonic regularization 



As expected, the propagators D^^ri Dy^y^ and Dy^„ derived in the previous section 
contain a pole corresponding to the Higgs boson. Apart from this, there is of course 
also a tachyonic pole in these expressions. Its Euclidian position, s — — A^, is given 
by the following equation: 

and differs from the position of the Landau pole A^, = jjie^^'^^/'^ by a correction of 
order / K\. The occurrence of this Euchdian pole creates problems because it spoils 
the causality of the theory. 

Let us analyse the way one performs a calculation in the 1/A^ expansion. To 
calculate, let's say, the Higgs propagator, one first performs a double expansion of 
the self-energy in the coupling constant A and in Then each coefficient in the 

expansion is given by a power series in the coupling constant A. For one given 
1/N order, due to the combinatorial structure, one is able to calculate explicitly 
all Feynman diagrams in the perturbative expansion, and to sum up the A series. 
The result then appears to contain an additional tachyonic pole. However, let us 
notice that the coefficient of the 1/A^ expansion which we calculate is not defined 
uniquely by the expansion in A which we calculate by Feynman diagrams. One still 
has the freedom to add a function in A whose perturbative expansion vanishes. Our 
treatment of the tachyonic pole is based on the observation that its residuum is indeed 
such a nonperturbative function of the coupling constant, of the type e^/^. One can 
add an opposite contribution, so that the tachyonic pole is cancelled, because this 
contribution vanishes exactly in all orders of perturbation theory. 

For this reason, we simply subtract the tachyon minimally at its pole from the 
leading order propagators. With this prescription, the calculation becomes well- 
defined, and one can represent the 0(1/A^) contribution in terms of a number of 
Feynman-like graphs. Such a tachyonic regularization can be carried out order by 
order in the 1/N expansion. 

Along the lines of the discussion above, we choose to regularize the propagators 
of cqns. 9 by subtracting the tachyon pole minimally. We introduce the following 
subtractions: 



^ ism\s) / 1 K \ 

n im^{s) / 1 \ 
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were 



1^ = (13) 

1 J ill 

is the residuum of the tachyonic pole in the aa propagator. 

It is worthwhile noting in this context that the tachyon not only has a mass of 
the order of the Landau pole mass. Also its residuum, i.e. its coupling to the other 
fields, is suppressed essentially by the same mass scale. Therefore the effects of the 
tachyon on the low energy physics are completely negligible as long as the mass of 
the Landau pole is not very low. 

Let us examine the tachyonic regularization of eqns. 12. When we use the tachy- 
onic regularized expressions instead of eqns. 9 for calculating higher order Green 
function, we essentially modify the Green function by a quantity proportional to k. 
As explained, this is a contribution which exactly vanishes in all orders of perturbation 
theory because k is a function of the type e^^^. 

Our tachyonic regularization can be seen as a different prescription for summing 
up the perturbative expansion of a given term in the l/A?" expansion. This prescription 
is such that it does not result in the presence of a tachyon in the theory's spectrum, 
which would violate causality. At the same time, the original information we started 
with, that is, the power series in A of the coefficient of the expansion, remains 
untouched. No doubt, this prescription is not unique. For a trivial theory, this am- 
biguity can be used for modeling the nondecoupling effects associated with unknown 
physics at a higher energy scale, very much in the same way a cutoff would do. One 
expects this uncertainty to become important numerically only when the energy of 
the process under consideration approaches the Landau scale. 

Finally, we would like to stress that we do not interpret the tachyonic regulariza- 
tion as a proof that the 1/N expansion is free of tachyons. It only says that within 
the usual derivation of the expansion the occurrence of tachyons is arbitrary. To 
establish conclusively the presence or the absence of tachyons, one needs a procedure 
for calculating the coefficients of the expansion without using Feynman diagrams 
at intermediary stages. Such a procedure is not available so far. Until substantial 
progress is made in this direction, we prefer a scheme which respects causality and 
allows one to perform higher order calculations consistently. 

4 Next-to-leading order calculation 

With the tachyonic regularized propagators defined at leading order in the previous 
section, it is now possible to give a diagrammatic representation of the next-to-leading 
order contributions to the self-energy functions. The graphs are given in fig. 2. Each 
graph is in fact an infinite sum of multiloop Feynman diagrams which all are of the 
same order in the expansion parameter 1/N. This is shown explicitly in fig. 3 for one 
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Figure 2: Infinite sums of multiloop Feynman diagrams which contribute in next-to- 
leading order in 1/N to the two-point functions of the 0{N) sigma model. The blob 
on propagators denotes the summed-up leading order propagators. Note that the tttt 
propagator at leading order in 1/N is a free propagator. One of the graphs above is 
shown in expanded form in fig. 3. 



graph. 

With the next-to-leading order self-energy functions defined in fig. 2, it is 
straightforward to calculate the propagators. In particular, the acr propagator, which 
is our main concern in this paper, is given by: 



m? (co'"'')( 



(14) 



where 



2 (corr) I 



s) 



m?(s) 



'Pis] 



(15) 



Because of the complicated form of the expressions, the evaluation of these graphs 
is highly involved and can only be performed numerically. In order to do the actual 
calculation new techniques were necessary. We developed techniques which extend 
the methods of ref. [|I^ which deals with massive three-loop diagrams. Some graphs 
are also ultraviolet divergent. It is rather complicated to disentangle the poles in 
e = n — 4 of the sets of multiloop diagrams of the type shown in fig. 2. Therefore in 



8 



3 

O 

1 Q' Q k 



^A{s)= { > - EE 

A:=0 1=0 



1 Q O I 

2 3 



Figure 3: Multiloop diagrams with three-loop topology which contribute to the XX 
propagator in next-to-leading order. 

order to perform renormalization one first determines the counterterms in terms of 
multiloop graphs by using the renormalization conditions discussed in the previous 
section. Then, without evaluating the counterterms explicitly, one adds them to the 
actual graphs which contribute to the self-energies a{s), -y{s) and S{s). In 

the next step, one identifies the divergences and sub divergencies of the graphs, and 
combines each of them with the appropriate terms of the counterterm contributions 
at the level of the integrands. In this way one renders each of the graphs ultraviolet 
finite. One can then apply the numerical methods of ref. [|l^] for calculating the 
ultraviolet finite graphs. 

At this point it becomes clear that identifying the divergences of the multiloop 
graphs, choosing an appropriate definition of the counterterms, imposing convenient 
renormalization conditions, and performing the actual numerical integration are in 
fact closely related issues. This is rather involved, and we will not enter into such 
details, which are beyond the scope of this article. We will discuss this in detail in 
a subsequent publication. Here we only give a sketch of the treatment for one of the 
graphs. 

Let us consider the graph Ai which is shown in fig. 3. This graph has the topology 
of a three-loop Feynman diagram, which is the most complicated topology among 
the graphs of fig. 2. To evaluate this graph, we treat it as if it were a three-loop 
diagram with two propagators replaced by some more complicated expressions, given 



by eqns. 12. We evaluate the resulting expression according to ref. |T^. This gives 
the following two-dimensional integral representation: 

/oo poo 
dp, / dpi^T, p^)D^^{P',)D^^{Pl)C\PlPlk') , (16) 
-oo JO 

where 



Po = pI-p' 

Pi = pi + 2poVi^+e-p' 

Pi = p2 + 2po/i + /i'-p' , (17) 
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and k = ^/s is the external momentum of the graph, fi is some arbitrary subtrac- 
tion scale for the subdivergencies of the graph, which will be used in the following. 
C{Pi,P2,Pi) = C{0,0,0;Pi,Pi,P^) is the three-point vertex diagram, for which 
we use an analytical relation given in ref. [0. That expression is suitable for our 
NLO calculation because it remains on the correct Riemann sheet for the range 
of parameters needed to evaluate the two-fold integral of eq. 16. 

The expression of eq. 16 cannot be evaluated numerically directly because it is 
ultraviolet divergent. First one has to combine it with the appropriate counterterms in 
order to subtract its divergencies and subdivergencies, thus rendering this expression 
finite. After some algebra, and after introducing the notation: 



1 



167r2p2 



log 



one finds the following subtracted expression: 



/oo poo , p -| 

dp, / dp (4vrp2) X {d^^{P^) [d^^{PI) - I^«(^o) i 



- 2D^^{P^,)D^^{Pl) [C{Pl Pi - C{Pl Pi p') + e] e 
+ 2D,,iPlD^,iP'^) [C{PI Pi e) - CiPl Pi /i^) + e] CiPl Pi /i^) 

+ D^xiPo)D^xiPi) \C{Pl Pi k') - C{Pl Pi p') + e 



(19) 



Here the subtraction scale of the subdivergencies, /i, can be set to zero without 
encountering any infrared difficulties. This choice is in fact more convenient for the 
actual calculation. The reader can easily convince himself that the expression above is 
ultraviolet convergent. The expression can be evaluated numerically fast and precisely 
by choosing an appropriate complex integration path, as explained in ref. [jT5| . 

The other graphs can be evaluated along the same lines. For the graphs A4 and A5 
we found it advantageous to rewrite the integrands in a special form, which results in 
a faster algorithm. These methods will be discussed in detail in future publications. 

At this point we would like to make some comments related to the renormaliza- 
tion procedure. First, we emphasize again that we use a renormalization scheme at 
intermediary stages, and a renormalization scale /i appears in relations like eq. 19. 
However, the final physical results are totally independent of this renormalization 
scheme, just like the leading order results. This is because the perturbation series of 
the 1/A^ coefficient is summed up exactly, to all orders. 

An interesting effect is that the wave function renormalization constants of the 
Goldstone and Higgs bosons become finite when calculated in the nonperturbative 
1/A^ expansion. This only shows up nontrivially in the NL order, as there are no 
contributions to the wave function renormalization constants in leading order. The 
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reader can easily convince himself that the only graphs potentially giving ultraviolet 
divergent contributions to the wave function renormalization constants are Bi and 
D. If one considers the Feynman diagrams which are contained in both Bi and D, 
one notices that they are all indeed ultraviolet divergent. However, after summing up 
all Feynman diagrams which compose these two graphs, the ultraviolet divergency of 
the wave functions disappears. One is left with very slowly (logarithmical) converging 
expressions. 

The finiteness of wave function renormalization constants is known to be a feature 
of the exact solution for interacting fields. In fact, one can prove they ought to 
be numbers with values between and 1. For the truncated perturbative solution 
this is notoriously not the case, wave function renormalization constants being even 
ultraviolet divergent in general. For the case of the NLO calculation described here, 
we were able to prove that the Higgs and Goldstone wave function renormalization 
constants are finite. However, we could not so far check that their values are between 
and 1. It is difficult to calculate their actual values numerically because their 
numerical convergence is only logarithmical and much slower than that of the graphs 
needed for physical calculations. 

The behaviour of higher orders in is an open question. Still, it is remarkable 
that the NLO of the expansion recovers such a property of the exact solution 
which is absent in perturbation theory. 



5 Nonperturbative Higgs lineshape 

The Higgs propagator which we calculated in the previous section is essentially a 
physical quantity. It can in principle be measured in a fermion scattering process, for 
instance the /i/i ti process, which is relevant for muon collider studies. The only 
contribution which one still has to add is the correction to the two Yukawa couplings 
of the Higgs to the fermions. As long as one only considers the contributions of 



leading order in the mass of the fermions, it is well-known [jT2|, |T3|, |1^ that these 
contributions are energy independent, and are simply given by the ratio of the wave 
function renormalization constants of the Higgs and the Goldstone bosons. Regarding 
the evaluation of this ratio, we already mentioned in the previous section that and 

are both finite at NLO in but evaluating them separately is very difficult 

because they converge very slowly. However, they have the same ultraviolet behaviour 
and for this reason their ratio is well-behaved numerically. In this context one can 
note that this ratio is ultraviolet finite even in perturbation theory. 

The results for the lineshapes are given in fig. 4. We plot the lineshape function 
/(^i) = \D„„{s)Zl^/Zl^\^, both at LO and at NLO in the expansion. We also 
compare with the known perturbative results up to NNLO [|T^. For the purpose 
of comparing the predictions of perturbation theory and of the expansion, we 
compare in fig. 4 Higgs lineshapes which have the same position of the peak. 

As one can see, there is a substantial discrepancy between the expansion in 
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Figure 4: The Higgs lineshape in fermion scattering. Comparison between the per- 
turhative lineshape (LO, NLO, NNLO) and the nonperturhative 1/N expansion (LO, 
NLO) for different values of the coupling. We compare lineshapes so that the position 
of the peak in two-loop perturbation theory and NLO 1/N expansion is the same. The 
lineshapes in NNLO perturbation theory and NLO 1/N are almost indistinguishable 
up to about 700 GeV. 
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leading order and perturbation theory even for low values of the coupling. We notice 
that our NLO result lies between the perturbative results and the lowest order in 
1/iV. The corrections are actually quite large compared to the lowest order in 
For low coupling we see that one is very close to perturbation theory now, which 
indicates that we should have a good approximation for the Higgs hue shape. For 
Higgs masses larger than about 850 GeV there is a fairly large difference between 
perturbation theory and the 1/N expansion. Perturbation theory is probably not 
reliable anymore, but things seem to be converging. 

For even larger values of the coupling the behaviour of the 1/N solution and of 
perturbation theory becomes different, and lineshapes cannot be compared at the 
same position of the peak any longer. This behaviour can be seen in fig. 5. Here 
we plot lineshapes in NLO expansion and in two-loop perturbation theory for 
a range of couplings. Since perturbation theory breaks down for large couplings, the 
perturbative plots for heavy Higgs bosons are given only for comparison purposes. 
While the position of the peak in perturbation theory increases continuously with the 
couphng, the 1/N solution shows a saturation around 930 GeV, after which only the 
decay width continues to increase. 

In order to study this in somewhat more detail, we plotted in fig. 6 two observable 
quantities which can be interpreted as an effective Higgs mass and width. Here, 
ideally one would hke to plot instead the real and the imaginary part of the pole of 
the Higgs propagator because they are universal for all Green functions and are not 
process dependent. However our numerical programs do not allow us at the moment 
to calculate the graphs on the second Ricmann sheet and solve the pole equation. 

Therefore we use a procedure to define an effective mass Mp^ak and width Fpeafc 
from the hneshape itself, which is more directly related to experiment. The mass 
Mpeak is defined as the location of the peak of the cross section. The width Fpgafe is 
taken from the peak height of the cross section as if the lineshape was determined by 
a resonant Higgs propagator only, giving rise to a Breit-Wigner shape: 

{S IVlpeak) J^V^peak^ peak 

Of course, the lineshape is not exactly a Breit-Wigner resonance. However, we 
found that this definition of the width is within a few percent of for instance the width 
at half-maximum. The description by Mpeak and Fpeafe gives therefore a satisfactory 
approximation of the leading features of the Higgs lineshape. It should be good 
enough to be used as a first approximation for phenomenological applications. 

Our results can strictly speaking not be used directly at the LHC. This is because 
one should also take higher order correction in the Higgs- Goldstone boson vertex 
into account. However, we expect that the correction to the width will be the dom- 
inant feature in practice. For values of the Higgs mass larger than about 900 GeV, 
the saturation effect sets in, and the uncertainty due to the discrepancy between the 
perturbative and the 1/N widths becomes important. The Higgs becomes wider, and 
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Figure 5: Higgs lineshapes in perturbation theory (two-loop) and in the nonperturba- 
tive 1/N expansion (NLO), where we marked the position of the peak. One can see 
the saturation of the position of the peak for large values of the coupling. 
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Figure 6: The relation between the peak variables Mpeak and Tpeak in perturbation 
theory and in the 1/N expansion. 
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the signal is washed out and lost in the background. This saturation effect can be 
quite important for LHC physics and the discovery limits there. 

We plot in fig. 6 the effective width Fpeafc as a function of Mpeak both in the 
expansion and in perturbation theory. One can see from fig. 6 that the expansion 
and perturbation theory appear to be converging towards a common relation between 
the Higgs width and mass. For a Higgs mass up to about 850 GeV, the agreement is 
remarkable. Afterwards the numerical values start to deviate. The nonperturbative 
width is larger than the width in perturbation theory, and saturation sets in. A 
maximum of the effective Higgs mass appears at around 930 GeV. 

For phenomenological purposes we give in the following a simple approximate 
formula to relate Mp^ak and Tp^ak- Expressing Mp^ak and Fpeafc in TeV, one can use 
the following relation: 

Ml^^ = 2.02582 ■ F,,,, - 0.812725 • Fj,,, 

-1.01729 ■F5,,, + 0.541511 -Fj,,, , (21) 

which is simply a fit of the NLO 1/A^ curve shown in fig. 6. Its precision is at 
the per mille level for the range of couplings shown in fig. 6. For comparison, the 
corresponding fit for the two-loop perturbative curve is the following (the exact two- 
loop lineshape is given in ref. p!6[]): 



M;,,, = 2.09069 -Fpeafc- 1.23887 -F^,,, 

peak + 0.115571 ■ J- pg^^ 



0.0233177 ■ + 0.115571 ■ . (22) 



6 Conclusions 

We treated an 0(A^)-symmetrical sigma model with spontaneous symmetry breaking 
in the nonperturbative 1/A^ expansion at next-to-leading order. 

In order to calculate corrections of higher order in 1/iV one needs a prescription for 
the treatment of the tachyons which appear in the propagators at leading order. We 
showed that in the usual derivation of the 1/A^ expansion the occurrence of tachyons 
is arbitrary, and it can be eliminated when one realizes that the coefficient of the 
1/A^ expansion is only determined up to a function of the coupling constant which 
vanishes identically in perturbation theory. Therefore we introduced a scheme of 
tachyonic regularization which can be performed consistently order by order in the 
1/A^ expansion. Our tachyonic regularization can be interpreted as a prescription for 
the summation of the coefficient of the 1/A^ expansion which leaves unchanged its 
perturbative expansion and preserves causality at the same time. 

The tachyonic regularization is not unique. This arbitrariness can be used for 
modelling nondecoupling effects of unknown physics at a higher energy scale. This 
arbitrariness is an interesting problem which deserves further investigation. 
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We then developed techniques to calculate infinite sets of multiloop Feynman 
diagrams which are needed for higher order calculations in 1/N. We applied these 
techniques to the Higgs sector of the standard model, = 4, and calculated the Higgs 
propagator and the lineshape of the Higgs resonance at muon colliders. The results 
are interesting, and confirm qualitatively some results obtained previously based on 
the leading order solution, such as the saturation of the Higgs mass. In NLO, it 
occurs at a value of about 930 GcV. Moreover, compared to the leading order, the 
NLO correction is rather substantial. It leads for low couplings to a remarkable 
agreement between perturbation theory and \/N expansion, thus solving a long- 
standing puzzle. This agreement proves that the nonperturbative 1/N expansion in 
higher orders, when combined with the tachyonic regularization, is a useful tool for 
treating strong couplings. It can give very precise results. 

Our NLO solution has certain properties associated with a nonperturbative solu- 
tion. For instance, it contains no residual renormalization scheme dependence. Of 
course, this is true for the leading order as well. An interesting property which only 
shows up at NLO is the finiteness of wave function renormalization constants, which 
also ought to be satisfied by the exact solution, and which does not hold in pertur- 
bation theory. 

Regarding Higgs searces at future colliders, it is interesting that the perturbative 
two-loop correction is substantial when compared to the one-loop and tree level, but 
turns out to be close to the nonperturbative solution which we derived up to about 
900 GeV. This reduces drastically the theoretical uncertainty of a heavy Higgs width. 

For larger couplings, our results show a deviation from perturbation theory. Es- 
sentially, the Higgs width is always larger than in perturbation theory. The Higgs 
mass is saturated, and its width can grow without its mass becoming larger at the 
same time. This effect must be taken into account in view of heavy Higgs searches at 
LHC. For phenomenological purposes we gave a simple approximate formula to relate 
the Higgs mass and width which is valid for large couplings as well. 
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